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IT) \ Abstract 

o 
o 



The high precision and large kinematic coverage of the data from the HERA-I 
^Sj ■ running period (1994-2000) have already allowed precise extractions of proton 

O . parton distribution functions. The HERA-II running program is now underway 

<D \ and is expected to provide a substantial increase in the luminosity collected at 

HERA. In this paper, a study is presented which investigates the potential im- 



pact of future data from HERA on the proton PDF uncertainties, within the 

currently planned running scenario. Next-to-leading order QCD predictions 

C\| ' for inclusive jet cross sections at the LHC centre-of-mass energy are presented 

^ . using the estimated PDFs. Finally, the effect of a possible future measurement 

| of the longitudinal structure function, F^, on the gluon distribution is investi- 

(N ' gated. 
On ' 
O 

■ 1 Introduction 

o 

Since the advent of HERA, much progress has been made in determining the Parton Distribution Func- 
tions (PDFs) of the proton. A knowledge of the PDFs is vital in order to make predictions for both 
^ ■ Standard Model and beyond the Standard Model processes at hadronic colliders. Furthermore, they must 

\ be known as precisely as possible in order to maximise the discovery potential for new physics at the 

LHC. HERA is now in its second stage of operation. With the measurements that can now be expected 
from HERA-II, our knowledge of PDFs should be further improved. In this paper, first studies of the 
potential impact on the PDF uncertainties of future measurements from HERA, are presented. 



2 HERA physics and kinematics 

Lepton-proton deep inelastic scattering (DIS) can proceed either via the neutral current (NC) interaction 
(through the exchange of a 7* or Z°), or via the charged current (CC) interaction (through the exchange 
of a W 1 * 1 ). The kinematics of lepton-proton DIS are described in terms of the Bjorken scaling variable, 
x, the negative invariant mass squared of the exchanged vector boson, Q 2 , and the fraction of energy 
transferred from the lepton to the hadron system, y. The three quantities are related by Q 2 = sxy, where 
s is the centre-of-mass energy squared. 

At leading order (LO) in the electroweak interaction, the double differential cross section for the 
NC DIS process is given in terms of proton structure functions, 

d2 ^dQ2 P) = %: [ y + F ^> Q 2 ) ~ Q 2 ) t Q 2 )] (i) 

where Y± = 1±(1— y) 2 . The structure functions are directly related to the PDFs and their Q 2 dependence 
is predicted by perturbative QCD. In particular, F2 and xF% depend directly on the quark distributions. 
F° r Q 2 ^5 1000 GeV 2 , F2 dominates the lepton-proton scattering cross section and for x < 10 -2 , F2 
itself is dominated by sea quarks while the Q 2 dependence is driven by gluon radiation. Therefore, HERA 
data provide vital information on the sea-quarks and gluon at low-x. At high Q 2 > M|, the contribution 
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Table 1 : The data-sets included in the ZEUS-JETS and HERA-II projected PDF fits. The first column lists the type of data 
and the second gives the kinematic coverage. The third column gives the integrated luminosities of the HERA-I measurements 
included in the ZEUS-JETS fit. The fourth column gives the luminosities assumed in the HERA-II projection. Note that the 
96-97 NC and the 94-97 CC measurements have not had their luminosity scaled for the HERA-II projection. 



from XF3 becomes increasingly significant and gives information on the valence quark distributions. The 
longitudinal structure function, Fl, is directly sensitive to the gluon, but is only important at high-y. 

At LO, the CC cross sections are given by, 

d 2 a cc (e+p) G 2 F ML r/ x2/ , Nl 

^dV^ = 2,x(Q 2 + M^) 2 X K fi + C_) + (1 " y) (d + S) ] 
d 2 a cc (e-p) G|Mi r/ . . \2 / t \i 

so that a measurement of the e+p and e~p cross sections provides information on the d- and u-valence 
quarks, respectively, thereby allowing the separation of flavour. 

The QCD scaling violations in the inclusive cross section data, namely the QCD Compton (7**7 — > 
gq) and boson-gluon-fusion (7*3 — > qq) processes, may also give rise to distinct jets in the final state. Jet 
cross sections therefore provide a direct constraint on the gluon through the boson-gluon-fusion process. 



3 PDF fits to HERA data 

PDFs are usually determined in global fits [1-3] made within the conventional DGLAP [4] formalism 
at next-to-leading order (NLO). Such fits use data from many different experiments, with the inclusive 
cross sections from DIS experiments providing the most important source of information. However, the 
high precision and wide kinematic coverage of the existing HERA data have already allowed precise 
extractions of the proton PDFs using only HERA data [5,6]. The high statistics HERA neutral current 
data determine the low-x sea and gluon distributions while information on the valence quark distributions 
is provided by the higher-Q 2 neutral and charged current data. The gluon PDF contributes only indirectly 
to the inclusive DIS cross section. However, it makes a direct contribution to jet cross sections. 

Recently, the ZEUS Collaboration have performed a combined NLO QCD fit to inclusive neutral 
and charged current DIS data [7-12] as well as high precision jet data in DIS [13] and 7p scattering [14]. 
This is called the ZEUS-JETS PDF fit [5]. The use of only HERA data eliminates the uncertainties from 
heavy-target corrections and removes the need for isospin symmetry assumptions. It also avoids the dif- 
ficulties that can sometimes arise from combining data-sets from several different experiments, thereby 
allowing a rigorous statistical treatment of the PDF uncertainties. Furthermore, PDF uncertainties from 
current global fits are, in general, limited by (irreducible) experimental systematics. In contrast, those 
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Fig. 1 : The optimised jet cross sections included in the HERA-II projected fit. The solid points show the simulated data 
generated using the NLO QCD programme of Frixione-Ridolfi, using the CTEQ5M1 proton and the AFG photon PDFs. The 
error bars show the statistical uncertainties, which correspond to 500 pb^ 1 of HERA data. Systematic uncertainties have been 
neglected. The dashed line shows the NLO QCD prediction using the ZEUS-S proton and AFG photon PDFs. The shaded band 
shows the contribution to the cross section uncertainty arising from the uncertainty in the gluon distribution in the proton. 
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Fig. 2: The fractional PDF uncertainties, as a function of x, for the u-valence, d-valence, sea-quark and gluon distributions at 
Q 2 = 1000 GeV 2 . The red shaded bands show the results of the ZEUS- JETS fit and the yellow shaded bands show the results 
of the HERA-II projected fit. 



from fits to HERA data alone, are largely limited by the statistical precision of existing measurements. 
Therefore, the impact of future data from HERA is likely to be most significant in fits to only HERA 
data. 



3.1 The ZEUS NLO QCD fit 

The ZEUS-JETS PDF fit has been used as the basis for all results shown in this paper. The most important 
details of the fit are summarised here. A full description may be found elsewhere [5]. The fit includes 
the full set of ZEUS inclusive neutral and charged current e ± p data from HERA-I (1994-2000), as well 
as two sets of high precision jet data in e + p DIS (Q 2 >> 1 GeV 2 ) and jp (Q 2 ~ 0) scattering. The 
inclusive data used in the fit, span the kinematic range 6.3 x 10~ 5 < x < 0.65 and 2.7 < Q 2 < 30000 
GeV 2 . 

The PDFs are obtained by solving the NLO DGLAP equations within the MS scheme. These 
equations yield the PDFs at all values of Q 2 provided they are input as functions of x at some starting 
scale Qq. The resulting PDFs are convoluted with coefficient functions to give predictions for structure 
functions and, hence, cross sections. In the ZEUS fit, the xu v (x) (u- valence), xd v (x) (d-valence), xS(x) 
(total sea-quark), xg{x) (gluon) and x{d(x) — u{x)) PDFs are parameterised at a starting scale of Qq = 7 
GeV 2 by the form, 

xf{x) = pix P2 {l - x) Ps P(x), (3) 

where P(x) = (1 + p±x). No advantage in the x 2 results from using more complex polynomial forms. 
The normalisation parameters, p\(u v ) and pi(d v ), are constrained by quark number sum rules while 
pi (g) is constrained by the momentum sum rule. Since there is no information to constrain any difference 
in the low-x behaviour of the u- and d-valence quarks, P2(u v ) has been set equal to P2{d v ). The data 
from HERA are currently less precise than the fixed target data in the high-x regime. Therefore, the high- 
x sea and gluon distributions are not well constrained in current fits to HERA data alone. To account 
for this, the sea shape has been restricted by setting Pi(S) = 0. The high-ir gluon shape is constrained 
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Fig. 3: NLO QCD inclusive jet cross section predictions at ^/i=14 TeV in three regions of pseudo-rapidity. The yellow and 
blue bands show the PDF uncertainties from the ZEUS-JETS and HERA-II projected fits, respectively. 



by the inclusion of HERA jet data. In fits to only HERA data, there is no information on the shape of 
d — u. Therefore, this distribution has its shape fixed consistent with Drell-Yan data and its normalisation 
set consistent with the size of the Gottfried sum rule violation. A suppression of the strange sea with 
respect to the non-strange sea of a factor of 2 at Qq is also imposed, consistent with neutrino induced 
dimuon data from CCFR. The value of the strong coupling has been fixed to ot s {Mi) = 0.1180. After all 
constraints, the ZEUS-JETS fit has 1 1 free parameters. Heavy quarks were treated in the variable flavour 
number scheme of Thorne & Roberts [15]. Full account was taken of correlated experimental systematic 
uncertainties, using the Offset Method [3, 16]. 

4 Results and discussion 

In this section, the results of two separate studies are presented. The first study provides an estimate 
of how well the PDF uncertainties may be known by the end of HERA-II, within the currently planned 
running scenario, while the second study investigates the impact of a future HERA measurement of Fl 
on the gluon distribution. All results presented, are based on the recent ZEUS-JETS PDF analysis [5], as 
described in the previous section. 

4.1 PDF uncertainty estimates for the end of HERA running 

The data from HERA-I are already very precise and cover a wide kinematic region. However, HERA- 
II is now running efficiently and is expected to provide a substantial increase in luminosity. Current 
estimates suggest that, by the end of HERA running (in mid-2007), an integrated luminosity of 700 
pb" 1 should be achievable. This will allow more precise measurements of cross sections that are curently 
statistically limited: in particular, the high-Q 2 NC and CC data, as well as high-Q 2 and/or high-i?T jet 
data. In addition to the simple increase in luminosity, recent studies [17] have shown that future jet cross 



section measurements, in kinematic regions optimised for sensitivity to PDFs, should have a significant 
impact on the gluon uncertainties. In this paper, the effect on the PDF uncertainties, of both the higher 
precision expected from HERA-II and the possibility of optimised jet cross section measurements, has 
been estimated in a new QCD fit. This fit will be referred to as the "HERA-II projection". 

In the HERA-II projected fit, the statistical uncertainties on the currently available HERA-I data 
have been reduced. For the high-Q 2 inclusive data, a total integrated luminosity of 700 pb" 1 was as- 
sumed, equally divided between e + and e~. For the jet data, an integrated luminosity of 500 pb~ x was 
assumed. The central values and systematic uncertainties were taken from the published data in each 
case. In addition to the assumed increase in precision of the measurements, a set of optimised jet cross 
sections were also included, for forward dijets in collisions, as defined in a recent study [17]. Since 
no real data are yet available, simulated points were generated using the NLO QCD program of Frixione- 
Ridolfi [18], using the CTEQ5M1 [19] proton and AFG [20] photon PDFs. The statistical uncertainties 
were taken to correspond to 500 pb -1 . For this study, systematic uncertainties on the optimised jet cross 
sections were ignored. The simulated optimised jet cross section points, compared to the predictions of 
NLO QCD using the ZEUS-S proton PDF [21], are shown in Fig.Q 

Table [^lists the data-sets included in the ZEUS-JETS and HERA-II projected fits. The luminosi- 
ties of the (real) HERA-I measurements and those assumed for the HERA-II projection are also given. 

The results are summarised in Fig. 13 which shows the fractional PDF uncertainties, for the u- and 
d- valence, sea-quark and gluon distributions, at Q 2 = 1000 GeV 2 . The yellow bands show the results 
of the ZEUS -JETS fit while the red bands show those for the HERA-II projection. Note that the same 
general features are observed for all values of Q 2 . In fits to only HERA data, the information on the 
valence quarks comes from the high-Q 2 NC and CC cross sections. The increased statistical precision 
of the high-Q 2 data, as assumed in the HERA-II projected fit, gives a significant improvement in the 
valence uncertainties over the whole range of x. For the sea quarks, a significant improvement in the 
uncertainties at high-x is also observed. In contrast, the low-x uncertainties are not visibly reduced. This 
is due to the fact that the data constraining the low-x region tends to be at lower-Q 2 , which are already 
systematically limited. This is also the reason why the low-x gluon uncertainties are not significantly 
reduced. However, the mid-to-high-x gluon, which is constrained by the jet data, is much improved in 
the HERA-II projected fit. Note that about half of the observed reduction in the gluon uncertainties is 
due to the inclusion of the simulated optimised jet cross sections. 

4.1.1 Inclusive jet cross sections at the LHC 

The improvement to the high-x partons, observed in the HERA-II projection compared to the ZEUS- 
JETS fit, will be particularly relevant for high-scale physics at the LHC. This is illustrated in Fig. |3j 
which shows NLO QCD predictions from the JETRAD [22] programme for inclusive jet production at 
y/s = 14 TeV. The results are shown for both the ZEUS-JETS and the HERA-II projected PDFs. The 
uncertainties on the cross sections, resulting from the PDFs, have been calculated using the LHAPDF 
interface [23]. For the ZEUS-JETS PDF, the uncertainty reaches ~ 50% at central pseudo-rapidities, 
for the highest jet transverse energies shown. The prediction using the HERA-II projected PDF shows a 
marked improvement at high jet tranverse energy. 

4.2 Impact of a future HERA measurement of Fl on the gluon PDF 

The longitudinal structure function, Fl, is directly related to the gluon density in the proton. In principle, 
Fl can be extracted by measuring the NC DIS cross section at fixed x and Q 2 , for different values of y 
(see Eqn. A precision measurement could be achieved by varying the centre-of-mass energy, since 
s = Q 2 /xy « AE e E p , where E e and E p are the electron and proton beam energies, respectively. 
Studies [24] have shown that this would be most efficiently achieved by changing the proton beam 
energy. However, such a measurement has not yet been performed at HERA. 

There are several reasons why a measurement of Fl at low-x could be important. The gluon 
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Fig. 4: The gluon PDFs, showing also the fractional uncertainty, for fits with and without inclusion of the simulated Fl data, 
for Q 2 — 1.5, 5, 10 and 20 GeV 2 . The red shaded bands show the results of the ZEUS-JETS fit and the yellow shaded band 
show the results of the ZEUS-JETS+F L fit. 

density is not well known at low-x and so different PDF parameterisations can give quite different pre- 
dictions for Fl at low-x. Therefore, a precise measurement of the longitudinal sturcture function could 
both pin down the gluon PDF and reduce its uncertainties. Furthermore, predictions of Fl also depend 
upon the nature of the underlying theory (e.g. order in QCD, resummed calculation etc). Therefore, a 
measurement of Fl could also help to discriminate between different theoretical models. 

4.2.1 Impact on the gluon PDF uncertainties 

The impact of a possible future HERA measurement of Fl on the gluon PDF uncertainties has been 
investigated, using a set of simulated Fl data-points [24]. The simulation was performed using the 
GRV94 [25] proton PDF for the central values, and assuming E e = 27.6 GeV and E p = 920, 575, 465 
and 400 GeV, with luminosities of 10, 5, 3 and 2 pb^ 1 , respectively. Assuming that the luminosity 
scales simply as E 2 this scenario would nominally cost 35 pb" 1 of luminosity under standard HERA 
conditions. However, this estimate takes no account of time taken for optimisation of the machine with 
each change in E p , which could be considerable. The systematic uncertainties on the simulated data- 
points were calculated assuming a ~ 2% precision on the inclusive NC cross section measurement. A 
more comprehensive description of the simulated data is given elsewhere [24]. 

The simulated data were included in the ZEUS-JETS fit. Figure|4]shows the gluon distribution and 
fractional uncertainties for fits with and without inclusion of the simulated Fl data. The results indicate 
that the gluon uncertainties are reduced at low-x, but the improvement is only significant at relatively 
low Q 2 < 20 GeV 2 . 

4.2.2 Discrimination between theoretical models 

In order to assess whether a HERA measurement of Fl could discriminate between theoretical models, 
two more sets of Fl data-points have been simulated [26], using different theoretical assumptions. The 
first of the two sets was generated using the MRSG95 [27] proton PDF, which has a large gluon density. 
The PDFs were then convoluted with the NNLO order coefficient functions, which are large and positive. 
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Fig. 5: The distribution of the longitudinal structure function Fl at Q 2 —5, 10 and 20 GeV 2 . The blue, red and green points 
show the simulated Fl data-points, respectively labelled maximum, middle and minimum in Table|5| The blue, red and green 
shaded bands show the NLO QCD predictions, in the case where the data-points of the corresponding colour have been included 
in the fit. For comparison, the yellow shaded band shows the prediction of the ZEUS-JETS fit. 
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Table 2: Summary of the PDFs used to generate the simulated Fl data-points. The extreme maximum Fl points were generated 
using the MRSG95 PDF, and convoluted with NNLO coefficient functions. The middle points were generated using the GRV94 
PDF, and the extreme minimum points were generated using the MRST2003 PDF, which has a negative gluon at low-x. 



This gives the "maximum" set of Fl data-points. In contrast, the second set has been generated using 
the MRST2003 [28] proton PDF, which has a negative gluon at low-x and low-Q 2 , thus providing a 
"minimum" set of Fl data. The original set of Fl points described in the previous subsection lies 
between these two extremes. The details of all three sets are summarised in Tab. 13 

Figure|5]shows the results of including, individually, each set of simulated Fl data into the ZEUS 
NLO QCD fit. The results show that the NLO fit is relatively stable to the inclusion of the extreme sets 
of data. This indicates that a measurement of Fl could discriminate between certain theoretical models. 
However, it should be noted that the maximum and minimum models studied here were chosen specifi- 
cally to give the widest possible variation in Fl- There are many other alternatives that would lie between 
these extremes and the ability of an Fl measurement to discriminate between them would depend both 
on the experimental precision of the measurement itself, as well as the theoretical uncertainties on the 
models being tested. 
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